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In this paper, a complete control architecture is proposed for the voltage conditioning stage of a hybrid 
power generation system composed of a Stirling engine coupled with a supercapacitor. Such a solar 
energy-based generation system aims at providing electricity to off-grid regions. The novelty of the pro¬ 
posed architecture is that it completely handles constraints on all the state variables of the electric stage 
while providing near to optimal performances in terms of settling time. The derivation of the control law 
enables a deep understanding of the main issues involved in the success of the closed-loop control. More¬ 
over, the resulting feedback laws are real-time compatible and are given in a complete explicit form. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Today, more than 1.4 billion people have no access to electric¬ 
ity, essentially in India, Indonesia, Bengladesh, Nigeria and Sub- 
Saharian Africa. Classical alternative fuel-based generators such 
as Diesel engines suffer from many drawbacks including their 
low reactivity to load changes, their high emissions and noise to 
mention but few issues [1 . 

Stirling engines are based on the Stirling thermodynamic cycle 
which describes a fluid cycling between a cold and a hot sources in 
order to transform heat energy into kinetic energy. This means that 
well isolated tanks containing water that is heated during the sun¬ 
ny hours can be viewed as green batteries that can provide kinetic 
energy to steer an alternator and produce electricity as soon as re¬ 
quired. In addition to this, the Stirling engines have the following 
advantages: low noise level, robustness, simplicity of the associ- 
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ated hydraulic circuit and a good thermal to electrical conversion 
efficiency. Then a micro solar power plant based on a Stirling en¬ 
gine seems to be a promising solution. 

However, one of the main issues of such a plant is its low re¬ 
sponse time when moving from an operating point to another 
when the supplied electrical load varies. An energy buffer which 
is a supercapacitor in the studied system is needed to palliate 
the load power demand in a fast manner while the Stirling engine 
adapts its produced power in a rather slower manner. Then an effi¬ 
cient control of the Stirling engine through its associated electrical 
part has a great influence on the dimensioning of the energy buffer 
that has to be as small as possible to reduce the cost of the plant. 

Within this framework, a control architecture for the voltage 
conditioning stage associated to a Stirling engine was developed 
and described in this paper. A preliminary solution to this problem 
has already been proposed in [2 where a one step unconstrained 
Model Predictive Control has been used. However, the approach 
proposed in [2] does not explicitly tackle the state constraints. In¬ 
stead an ad hoc filter has been added in order to enlarge the do¬ 
main of applicability of the proposed solution. 
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Many hybrid power generation systems have been studied in 
the literature either for grid or off-grid applications. They combine 
different renewable energy sources essentially wind, solar photo¬ 
voltaic, hydraulic and fuel cells. Different power electronics archi¬ 
tectures associated to those hybrid power systems have been 
proposed and their corresponding control strategies developed. 
In [3 , a variable speed wind turbine associated to a simple power 
electronics have been used to supply AC loads, an MPPT (Maximum 
Power Point Tracking) algorithm was used to extract maximum 
power from the wind using a DC/DC boost converter while an in¬ 
verter was used to create the AC grid and to regulate the DC bus, 
however because of the random nature of the wind, the AC voltage 
was perturbed. In [4], a similar system combined to a Diesel Engine 
and a battery bank was studied. The batteries act as a buffer source 
(through a buck-boost converter) to regulate the DC bus absorbing 
power in case of high wind and discharging in case of scarce wind. 
A linearized models were used to control the system at some oper¬ 
ating points. However for many applications where the loads and 
power sources can vary significantly, such controllers are not suit¬ 
able. In [5 , a similar system provided by dump loads and its asso¬ 
ciated switch was studied and PI based controllers were proposed. 
In [6] a hybrid wind/photovoltaic/fuel cell power system for grid- 
connected applications was proposed. An MPPT techniques were 
used to extract maximum power from the wind and PV panels 
while the fuel cells system suppress the grid power fluctuations. 
The objectives of such hybrid systems is to act as a battery-less 
uninteruptible power source (UPS).Valenciaga et al. ( 7-9 ) pro¬ 
posed an interesting technique to control a hybrid wind/photovol¬ 
taic (and even extra power sources) with wind as a primary source 
and a battery bank on the DC bus as an energy storage device. The 
control design proposed by Valenciaga et al. are based on passivity/ 
sliding mode techniques (and a second order sliding mode tech¬ 
nique in [9]) for different wind regimes. Supervisory controllers 
for such systems were proposed in [10-12]. All the techniques pre¬ 
viously cited do not handle explicitly the constraints of the sys¬ 
tem’s state variables and controllers which can lead to some 
undesirable behaviors. 

In the present paper, a complete state feedback solution is pro¬ 
posed that entirely solves the constrained problem while providing 
a deep understanding of the interaction between the different sub¬ 
systems composing the whole voltage conditioning stage. This is 
done by decomposing the whole system into subsystems for which 
the vector fields can be analyzed graphically. A nice feature in the 
proposed solution is that it suggests a general control design meth¬ 
odology for a family of power circuits despite the fact that its suc¬ 
cess may depend on the quantitative design of the circuit 
components. The proposed design underlines the conditions under 
which the framework is successful and can therefore be also used 
to guide the design step. 


The paper is organized as follows: First of all, the solar micro 
plant is described in Section 2. The mathemetical equations 
describing the system’s dynamics which are used in the design of 
the control law are described in Section 3. Some preliminary re¬ 
sults that are used in the remainder of the paper are derived in Sec¬ 
tion 4. The explicit expression of the feedback law is described in 
Section 5 and simulation-based validation is proposed in Section 
6 through a wide variety of scenarios in order to strengthen the 
successful constraint handling capability of the proposed feedback. 
Finally, Section 7 concludes the paper by summarizing its contribu¬ 
tions and by giving hints for future investigations. 

2. Description of the solar micro plant 

Fig. 1 shows an overview of the solar micro plant including the 
Stirling engine linked to a synchronous machine together with its 
voltage conditioning stage that consists of a diode bridge, a full 
bridge DC/DC converter, a DC bus, a single phase rectifier. A sup¬ 
ercapacitor is also present and linked to the DC bus through a bidi¬ 
rectional DC/DC converter. 

More precisely, the Stirling engine is a thermodynamic machine 
that transforms heat (of the hot fluid flowing inside one side of the 
Stirling engine internal heat exchanger) into a mechanical work by 
performing a Stirling thermodynamic cycle at each shaft rotation 
[13]. The produced mechanical torque is used to drive the synchro¬ 
nous machine that reacts by producing an electromagnetic torque 
which together with the motor torque impose the rotational speed 
of the shaft. The variable speed of the engine induces a variable 
amplitude and frequency voltage at the output of the PMSG which 
is conditioned through the mentioned power electronics stage be¬ 
fore feeding the load with a 230 V and 50 Hz voltage. The process 
of the voltage conditioning is as follow: 

The variable amplitude and frequency voltage at the output of 
the PMSG is rectified through the diode bridge to produce a vari¬ 
able DC voltage which is conditioned through a DC/DC full bridge 
converter that is connected to a DC bus whose voltage V bus needs 
to be regulated at the appropriate value of 50 V which is compati¬ 
ble with the good functioning of the inverter that supplies the load 
with the AC voltage. 

Note however that in order to dynamically meet the load de¬ 
mand, the Stirling engine alone is not sufficient due to its high 
mechanical and thermal inertia. That explains the use of the sup¬ 
ercapacitor (see Fig. 1) as a fast energy buffer during highly dy¬ 
namic transients to tightly regulate the DC bus voltage V bus at 
50 V through appropriate management of the control u 2 (duty ratio 
of the bidirectional DC/DC converter connected to the supercapac¬ 
itor). The control iq (duty ratio of the DC/DC full bridge converter) 
has then to be used to control the output current of the DC/DC full 
bridge converter namely I Lfb at some desired reference value 
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Fig. 1 . Overview of the solar micro plant including the Stirling engine together with its voltage conditioning stage. 
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computed so as to meet the load power demand and to enable the 
charge/discharge of the supercapacitor in order to recover its de¬ 
sired voltage V r s f = 120 V. 

Note that the engine torque controller depicted at the left edge 
of Fig. 1 is not studied in this contribution. Indeed, this controller 
uses the information on the power demand, the input tempera¬ 
tures of the hot and cold fluids flowing through the internal heat 
exchangers of the Stirling engine (see [13 ) in order to compute 
the mass flow rates set-points leading to an optimal thermoelectric 
efficiency [14]. This loop, while important for the global efficiency 
of the systems, obviously works in a separated time scale when 
compared to the characteristic of the voltage conditioning stage. 
This makes relevant a separated control design of the latter that 
is based on the assumption of a constant torque at the input of 
the PMSG block. 


3. Mathematical model and problem statement 

The Stirling engine studied in this paper is a prototype available 
within Schneider Electric. It is a 3 kW e/ec (maximum electrical out¬ 
put at 300 °C for the hot fluid temperature) alpha type Stirling en¬ 
gine (see [13]). From steady state experiments performed on this 
engine, an expression of the motor torque was fitted to the exper¬ 
imental measurements using the Eureqa-Formulize® equation fit¬ 
ter ( 15,16 ) leading to the following result: 

Tmot = a (m h ,m c , f”, T?) ■ Q(t) + j?(ih h , m c , T ”, f c n ) (1) 

where T mot is the Stirling engine output torque, Q(t) is the shaft 
speed (in rad/s), m h and m c are the mass flow rates of the hot and 
cold fluids at the inlet of the Stirling engine, Tj] 1 and T" are their 
respective temperatures. a(.) and /?(.) are the functions fitted with 
the experimental measurements and given by: 

a(m h , m c , T]? jf)= 0.2498 + 0.004109 • 1™ + 3.13 x 10” 5 • T% ■ T[ n • m h 

-1.81 2 • m h - 1.655 x 10” 5 ■ if • f c n - 0.007562 ■ f c n m h 
=0.1 T" 1 -30.03 

Fig. 2 compares the Stirling engine torque computed using Eq. 
(1) with the torque derived from the experimental data for a set 
of 160 operating points defined by m h , rh c , Tjj 1 and T™. It can be 
seen from Fig. 2 that the curve corresponding to Eq. (1) fits most 
of the experimental data points. 


In a recent contribution [2], it has been shown that the dynamic 
system lying to the right of the Stirling engine (see Fig. 1) can be 
described by the following set of ordinary differential equations: 


X\ = — cqxi — U3X2 + cl 2 

(2a) 

x 2 = -a 4 x 2 + a 6 x 1 - a 7 x 3 

(2b) 

x 3 = a 8 (x 2 - kx 4 U\) 

(2c) 

x 4 = a 9 (-x 5 +toc 3 Ui) 

(2d) 

x 5 = a w (x 4 +x 6 -- — L —) 

(2e) 

V ^1 vnv ' ^5 J 


x 6 = a n (-x 5 +x 7 u 2 ) 

(2f) 

X 7 = —0,\ 2 X§U 2 

(2g) 

where the state components x,, i e {1,.. 

,., 7} are associated to the 

physical variables depicted in Fig. 1 as follows: Xi = £2, x 2 = / re d, 


*3 = V red , x 4 = hfb, *5 = V bus , x 6 = I Lbb and x 7 = V sc . The coefficients 
at, i e { 1 ,..., 12 } are constant parameters depending on the com¬ 
ponent values according to the expressions given in [2]. rj inv in¬ 
volved in (2e) represents an efficiency factor. 

Note that the coefficients <q and a 2 depend on the mechanical 
torque of the Stirling engine defined by Eq. (1) and on the equiva¬ 
lent moment of inertia of the rotating shaft (see [2]). According to 
the comment made at the end of the preceding section, this coeffi¬ 
cients are supposed here to be constant as the torque variation is 
assumed to be very slow when compared to the electric stage char¬ 
acteristic times. 

Note also that the Eqs. (2a), (2g), (2c)-(2f) and (2g) describe the 
mean behavior so that iq and u 2 may take all values in [0,1]. The 
resulting control values are then implemented through standard 
Pulse Width Modulation (PWM) techniques. 

Note finally that the model used for the supercapacitor is sim¬ 
ply that of a big capacity whose value being equal to 5F neglecting 
the equivalent resistance in series since it is generally very small 
(as it is the case for the real supercapacitor available within Schnei¬ 
der Electric testbed). 

As far as numerical values are concerned, the coefficients 
a i: i e { 1 ,..., 12} for the studied system and for some operating 
point are given by: <q = -0.183, a 2 = 558.11, a 3 = 118.4453, 
a 4 = 9615.4, a 5 = 1.3712, a 6 = 5101.1, a 7 =641.02, a 8 = 425.53, 
a 9 = 6250, aio = 7.34, a u = 4484.3, a u = 0.2, rj inv = 0.95 and 
k= 1. 

Rephrasing the control objective in terms of the new state vari¬ 
ables leads to a control objective aiming at regulating x 5 and x 7 
around their desired values which are x r s = 50 V and x r 7 = 120 V 



Fig. 2. Comparison between the fitted Stirling engine torque equation and the torque derived from the experimental data, each point in the figure corresponds to an operating 
condition defined by m h , rh c , T™ and T™. 
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respectively. Moreover, the validity of the model together with the 
electrical requirement imposes the following set of constraints: 

X/G i e {2,3,4,6} (3) 

The constraints on the regulated variables x 5 and x 7 will be natu¬ 
rally respected by the success of the control loop while the con¬ 
trolled excursion on Xi is implicitly imposed through the 
constraints on both x 2 and x 3 as it can be inferred from Eqs. (2a) 

and (2b). 

4. Preliminary results 

In this section some preliminary results that are used in the se¬ 
quel are stated and, when needed, proved. Let us first begin by the 
following assumptions that are checked to be satisfied for the stud¬ 
ied system: 

Assumption 1. The linear system in the state (xi,x 2 ) defined by 
( 2 a), ( 2 b) (steered by the exogenous signal x 3 ) shows two 
separated negative real eigenvalues 2 / and 2 S such that I /L/ |>| I- 

For the studied system, this assumption is satisfied with 
| 2/ |« 9552 Hz >| |« 63.4 Hz. This feature is likely to be satis¬ 

fied on any testbed as it reflects the co-existence of mechanical 
(xi) and electrical (x 2 ) variables which generally evolve among sep¬ 
arated time scales. 

Assumption 2. The bounds x% in and x^ ax invoked in (3) are 
compatible with the constraints on x 2 , namely the steady values of 
x 2 belong to [x^ in ,x^ ax ] for all constant value of x 3 e [x™ n ,x^ ax ]. 
Moreover, x™ m > x r 5 . 

This assumption can be checked on Fig. 3 (see also Assumption 
3 below) where all possible steady states of (2a)-(2c) and (2d) in 
which x 5 = x^ is used, are depicted as functions of tq. It can be eas¬ 
ily checked that a wide variation range of x 3 corresponds to a small 
variation range of x 2 . In particular, Assumption 2 is clearly satisfied 
as soon as the interval [x]f n ,x™ ax ] contains the interval 
[4.736,4.919]. Note that the need for the condition x™ in > x r 5 can 
be easily understood by observing Eq. (2d) which indicates that 
the controllability of x 4 using iq e [ 0 , 1 ] (when x 5 ^x^) requires 
the inequality x 3 > x^ to be satisfied. 

The first result states that it is possible to monitor the fulfill¬ 
ment of the constraints on the rectified current x 2 = I red by 

x^(u\) (rad/s) 



*?(*!) (V) 



monitoring the excursion of the rectified voltage x 3 = V red . More 
precisely, the following lemma is derived: 

Lemma 1. There are easily computable maps x^Xi ,x 2 ) and x 3 (x\ , x 2 ) 
such that all trajectories of (2a), (2b) starting from an admissible 
initial conditions (regarding the constraints on x 2 ) and steered by a 
profile x 3 (-) meeting the following condition: 

X 3 (t) G [X 3 (X 1 (t),X2(t)),X3(X 1 (t),X 2 (t))] (4) 

satisfy the saturation constraints on x 2 . 

Proof. See appendix A where explicit expressions (A.2), (A.3) are 
given for x 3 (xi, x 2 ) and x 3 (xi, x 2 ). □ 

The next assumption concerns the dynamic of the subsystem 
described by Eqs. (2a)-(2c) and (2d) in which the state is given 
by z := (xi,x 2 ,x 3 ,x 4 ) r e U 4 while x 5 is assumed to be equal to its 
reference value, namely: x 5 = X 5 = 50V. Such a consideration is rel¬ 
evant as the tight regulation of x 5 around x^ will be achieved using 
the control u 2 as it is shown later. Obviously, under these condi¬ 
tions, the system ( 2 a)-( 2 c) and ( 2 d) can be written in the following 
condensed form: 

z=4(u,)z + bQ j (5) 

Assumption 3. For each positive value iq e [0,1], the matrix A(iq) 
is Hurwitz. Consequently, for each such constant value of iq, the 
system (5) shows steady states that are denoted hereafter by 
xf(iq), i e {1, • ■ • ,4}. 

The evolution of the steady state values xf(iq), i e {1 ,..., 4} 
are depicted in Fig. 3. This figure shows a very interesting feature 
regarding the evolution of xf(iq). Indeed, it is remarkable that a 
wide range of values of x 4 can be asymptotically produced (which 
is crucial in the sequel when it comes to regulate the super-capac¬ 
itor voltage x 7 = V sc ) for a rather very restricted range of values of 
x 2 (between 4.736 A and 4.919 A). In order to refer to this feature in 
the forthcoming argumentation, the following definition is used: 

Definition 1. The range of variations of steady values of x 2 is 
denoted hereafter by s 2 , namely 

c 2 := maxXo f (tq)- minx^(iq) ( 6 ) 

UlG[0,l] u lG [ 0,1] 


zf(ui) (A) 



xfM ( A ) 



Fig. 3. Illustration of Assumption 3: Evolution of the steady state values xf (ui), i e (1,... ,4} defined in Assumption 3 as a function of U\ for the studied system. 
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while the central value of xj(-) is denoted hereafter by xf c . Conse¬ 
quently, X^(iq) el 2 :=xf c + [-f ,+f]. 

Note that the parameters £ 2 and xf c are determined by the 
dimensioning of the system and does not depends on any control 
design. 

Fig. 3 shows that for the studied system, x s 2 c « 4.82 while 
£ 2 ^ 0.12 which corresponds to a relative variation of less than 
2.5%. 

Another feature that can be observed in Fig. 3 is that all the 
maps are monotonic so that the following notation can be used: 

xf(xf)X(*f) (U) c {1,...,4} 

to denote the stationary values of x, and iq that correspond to a gi¬ 
ven stationary value Xj of the jth component of the state vector z. 
This is extensively used in the sequel. 

Based on the observation made earlier regarding the excursion 
of the stationary values of x 2 , it can be argued that if the conclusion 
of Lemma 1 is used to define bounds on the variations of x 3 such 
that x 2 remains in a tight interval around x s 2 t,c (see Definition 1), 
then it becomes relevant to analyze the behavior of the dynamic 
system ( 2 c) and ( 2 d) in which the state is (x 3 ,x 4 ) while x 2 and x 5 
are considered to be constant. Such an analysis can then be per¬ 
formed by visualizing the vector fields in the plan (x 3 ,x 4 ). When 
doing this analysis using the numerical values mentioned at the 
end of Section 3, it comes that the following assumption is 
satisfied: 

Assumption 4. For any admissible desired pair (x^xj (X3)), the 
picture of the vector fields of the dynamic system ( 2 c)-( 2 d) in 
which x 2 = x 2 (x^) el 2 and x 5 = x r 5 is the one depicted in Fig. 4 
below, more precisely: 


1. The horizontal lines AB and DE correspond to the lines x 4 = x% in 
and x 4 = xj 11ax respectively. 

2. C is the desired steady point (x^xj (x^)) while u\ = uf (x^) is the 
corresponding steady control. 

3. DC is the trajectory under the control iq = 0 that starts at D and 
ends at C (this defines D). 

4. BC is the trajectory under the control 1/1 = 1 that starts at B and 
ends at C (this defines B). 

5. The disposition of the vector fields on the line AB and DE are 
such that the control iq (x) given by: 


Ui (x) := 


i^[X5+P 4 (X% ax -X 4 )} 

^[x 5 +P 4 (xf n -x 4 y 


on DE 
on AB 



steers the state towards D (resp. B) when x 4 = x% ax (resp. x% in ) while 
keeping x 4 constant. 

6 . The line ABCDE is attractive under the control iq = 1 for any 
starting point that lies to the right of the curve BCD (this region 
is denoted by 7Z 34 on Fig. 4). 

7. The line ABCDE is attractive under the control iq = 0 for any 
starting point that lies to the left of the curve BCD (this region 
is denoted by C 34 on Fig. 4) 


Note that the first 4 items of Assumption 4 are only definitions. 
The item (5) can be easily checked by investigating the disposition 
of the vector fields on the horizontal lines AB and DE. It is also pos¬ 
sible to formulate algebraic conditions [see B.l, B.2, B.3 and B.4] to 
check items (5), ( 6 ) and (7) above as it is shown in appendix B. 

Although the conditions B.l, B.2, B.3 and B.4 might seem tricky 
to check, it is likely that the topology of of the vector fields de¬ 
picted in Fig. 4 holds for a wide class of such voltage conditioning 
systems. The above conditions are only given for completeness. 

Note also that the disposition of the vector fields obviously sug¬ 
gests that a sliding modes controller can be defined in which the 


curves ABCDE is the sliding surface while the equivalent control 
on the sliding surface steers the system towards the desired point 
C. Flere we prefer to use continuous control because the Eqs. (2a)- 
( 2 c) and ( 2 d) represent the mean dynamic and in order to derive an 
easy to use closed-loop systems avoiding the numerical difficulties 
associated to the sliding modes controllers simulation when long 
tuning and assessing scenarios are required as it is the case here 
(see later). 

Rather than using the sliding modes formulation, a 2D interpo¬ 
lation techniques is used in the (x 3 ,x 4 ) plan. Namely, an interpo¬ 
lated expression of the feedback is derived in which, the 
expressions 1/1 = 1 or tq = 0 are obtained when the current state 
(x 3 ,x 4 ) is far from the curve ABCDE while iq takes the appropriate 
values when (x 3 ,x 4 ) lies on the curve ABCDE. The appropriate value 
is itself obtained by appropriate interpolation between the values 
depicted on Fig. 4 according to the position of the current point 
(x 3 ,x 4 ) on the curve ABCDE. 

In order to properly define the 2D interpolated expression, the 
following auxiliary function is needed: 

(p(r,v u v 2 ) := V\ + Sat(l - e^ r ) • (v 2 - v x ) (8) 

which simply takes the value V\ for r = 0 and v 2 when r > 0 with a 
transition rate defined by the positive constant L 

Based on the above discussion, the following explicit feedback 
can be defined for a given target point (x^xj(X 3 )): 

Definition 2. Let a desired point (x^x^x^)) be defined. Given a 
current position p := (x 3 ,x 4 ) such that x 4 e [xj^xj 0 *], compute the 
following quantities: 


• Compute u^ E (x) and iq B (x) according to: 

Uf (X) := <P , 0, a, (x) 

uf (x) := cp ( X3 ^ (B) X3 ,T Ul (X)) 

where tq(x) is given by (7). 

• Compute uf D (x) according to: 


/ / v max v 

X -X 4 


Ii^ CD (x) := < 


(p l 4 

<p 


X4-X 


min 


0 ,uj) 

if x 4 > x; 

i,ui) 

if x 4 % x; 


(9) 

( 10 ) 


where u\ := uf (x£). 

• Compute the abscissa x 3 of the horizontal projection of p on 
BCD. 



Fig. 4. Typical disposition of the vector fields of the dynamic system (lc) and (Id) 
when x 2 = x s { (x * 3 ) and x 5 = x r 5 are assumed to be constant. 
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Compute the vertical interpolation iif(x) according to: 


/ / v max v 

x. -x 4 


u"(x) ■- < 


<P 


(P 


,uf(x),l 


X4-X'. 


,uf(x),0 


if x 3 > x 3 
if x 3 ^ x 3 


*4" + *7 < - * 06) 

This obviously introduces constraints on the possible load power 
levels that are compatible with a sustainable operation. This is rig¬ 
orously stated in the following assumption: 


• Compute U\ (x | x^) by horizontal interpolation between uf(x) 
and u^ CD (x) according to: 

u,(x|x;) : =<?( |x j7* 3 ' ,urw,ur(x)) (ii) 

Note that the resulting feedback control iq(x | X 3 ) defined by 
( 11 ) is continuous with respect to x for a given target 
(X3,xj) = (x^xj(X3)). Indeed, for a given value of x 4 , continuity 
w.r.t x 3 is guaranteed by the continuity w.r.t x 3 of the expressions 
(9), (10) and (11) while the continuity w.r.t x 4 for a given x 3 is guar¬ 
anteed by the continuity w.r.t x 4 of the expressions of uf CD (x) and 
uf (x). 

Note that the success of the state feedback iq(x | X 3 ) in steering 
the state (x 3 ,x 4 ) to the desired (x^xj) is conditioned by the fact 
that the excursion of x 2 is small. Indeed, under this assumption, 
the picture depicted in Fig. 4 remains relevant while x 2 is slightly 
moving. Note however that for this excursion to be small, the fol¬ 
lowing two requirements are needed: 

• The range of variation s 2 of the steady values of x 2 (invoked in 

Definition 1) is small. 

• The variations of x 3 remains inside the moving interval 
te(Xi,x 2 ),x 3 (Xi,x 2 )] which is computed online based on the 
admissible range [x™’ n ,x™ ax ] [see (A.2), (A.3) in the appendix]. 

The first requirement mentioned above depends on the choice 
of the system’s components and has to be satisfied by construction. 
The second requirement can be satisfied by an appropriate man¬ 
agement of the desired value X 3 . This is explained in the following 
discussion. 

First of all, the following lemma is needed: 


Assumption 5. There is a positive real S 7 such that all possible 
values P L of the load power demand satisfy the inequalities (16). 

Before the control law u 2 = K 2 (x \ P L ) can be completely defined, 
the following lemma is needed: 


Lemma 3. The state feedback law given by: 



—— ■ (an • *5 - a 10 • (x 5 - x^) + x r 6 ef - p 6 • (x 6 

On • X 7 

-*?)) 


(17) 


where: 



Pl 

dinv ' ^5 


x 4 + ^ (x r 5 
Cl 10 




in which p 5 and p 6 are positive design parameters, locally asymp¬ 
totically steers the pair (x 5 (t),x 6 (t)) to the reference (x£,Xg(t)), 
namely: 



(19) 


Proof. See Appendix C. □ 

Now we have all we need to define the reference value X 3 that is 
used in Definition 2 to compute the state feedback i/i (x | X 3 ). This is 
defined to be the closest value to x s f (x 4 ) [see (14)] that is compat¬ 
ible with the requirement on the excursion of x 2 . More precisely: 

* 3 (*, p i) : = ar S min | z 3 -xf^(x, P L )) | (20) 

Z 3 G[X 3 (X),X 3 (X)] 


Lemma 2. Any appropriate control that asymptotically steers x 5 to its 
reference value x r 5 leads to the following property: 

lim[x 6 (f) -Xg(t)] = 0 (12) 

t—> oo 

where Xg(t) is given by: 

x£(f)=-A--x 4 (t) (13) 

w in v ' a 5 

where P L is the load power demand. 


Proof. This directly follows from Eq. (2e) in which x 5 is replaced by 
x^ while x 5 = 0 is used. □ 


It comes from Lemma 2 together with Eq. (2g) that it is possible 
to regulate the supercapacitor voltage x 7 by choosing the following 
reference value for x 4 : 


x 4 (x,Pl) 



+ <5 7 • Sign [xj) 



S 7 > 0 



since according to (12), (13) this asymptotically infers the following 
dynamics to x 7 : 

x 7 = ai 2 -u 2 -S 7 ■ Sign(x r 7 - x 7 ) (15) 


which obviously stabilizes x 7 around x r 7 . 

Note however that for (15) to be compatible with the con¬ 
straints on x 4 , the following conditions are needed: 


where x r 4 (x,P L ) is given by (14) depending on the current value of 
the supercapacitor voltage x 7 and the present load demand P L . Note 
that by doing so, the rate with which the supercapacitor voltage is 
regulated around its reference value is moderated by the dynamic 
of the mechanical part of the system which is directly linked to 
the Stirling engine performance. This is done through the definition 
of the admissible interval [x 3 (x 1 (t),x 2 (t)),x 3 (xi(t),x 2 (t))]. 

Remark 1. Note that in [2], it has been observed that in order to 
avoid negative excursions on some of the system’s variables, the 
variation of the set-point x r 4 needs to be filtered or slowed down. 
However, contrary to [2] where this feature has been observed , the 
discussion above gives a deeper interpretation and a more 
systematic, efficient and dynamic way to address this limitation. 
Indeed, while the filter introduced in [2] operates continuously 
slowing down the system unconditionally, the values that result 
from (20) slow down the system only when necessary. This is a 
general fact making explicit handling of the constraint preferable 
to worst-case-based unconstrained design. 

The preceding discussion completely defines the feedback law 

U, = K x (X,P l ) := U,(X I X* 3 (X,P l )) 

The remaining part of this section is dedicated to the definition of 
the feedback law u 2 =I< 2 (x,P L ). Fortunately, the same graphical 
analysis can be done in the 2D plan (x 5 ,x 6 ) for a given value of x 4 
[see Eqs. 2e,2f]. This analysis is relevant under '20) since as it is ex¬ 
plained earlier, Eq. (20) leads to rather slow variations of (x 3 ,x 4 ) 
inducing hence two separated time scales, a slow one for 
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Xi, i e { 1 , 2 ,3,4,7} on one hand and a fast one forx 5 and x 6 on the 
other hand. 

When considering the vector fields of (2e)-(2f) for a given value 
of x 4f it can be shown that the following assumption holds for the 
studied system: 


Assumption 6. For a given admissible x 4 and a supercapacitor 
voltage *7 >*5’ the disposition of the vector fields around the 
desired pair (x r 5 ,x r 6 ) where x’ 6 is defined by (13) is the one given in 
Fig. 5 below: more precisely: 


1. The horizontal lines AB and DE correspond to the lines x 6 = xg in 
and x 6 = x™ ax respectively. 

2. C is the desired steady point (xj^Xg) while u\ is the locally sta¬ 
bilizing feedback invoked in Lemma 3: 

3. DC is the trajectory under the control u 2 — 0 that starts at D and 
ends at C (this defines D). 

4. BC is the trajectory under the control u 2 = 1 that starts at B and 
ends at C (this defines B). 

5. The disposition of the vector fields on the line AB and DE are 
such that the controls u 2 given by: 


U 2 (x) := 


*7 

*7 


X5+fi(xr-*e) 

X 5+^0C"-*6) 


on DE 
on AB 



steers the state towards D (resp. B) when x 6 = x^ ax (resp. when 
x 6 = Xg !n ) while keeping x 6 constant. 

6 . The line ABCDE is attractive under the control u 2 = 1 for any 
starting point that lies to the left of the curve BCD (this region 
is denoted by C 56 on Fig. 5). 

7. The line ABCDE is attractive under the control u 2 = 0 for any 
starting point that lies to the right of the curve BCD (this region 
is denoted by 7 Z 5 g on Fig. 5). 


Note that the same comments that followed Assumption 4 can 
be reproduced here regarding the items invoked in Assumption 6 . 
In particular, there is here also the possibility to check them 
through explicit conditions similar to (B.1)-(B.3) and (B.4). 

Based on the observation of Fig. 5, an interpolation-based state 
feedback I< 2 (x | P L ) can be defined in a similar way to the one al¬ 
ready used in Definition 2 to define I< 1 (x | Pi). 

Definition 3. Let be given a load power demand P L . Consider the 
corresponding desired pair (x^Xg) using x ] 6 defined by (13). 
Compute the following quantities based on the current point 

p = (x 5 ,x 6 ) r : 


Xr 



Compute u% E (x) and irJ B (x) according to: 


AB, 


uf(x) -=(p( 


x 5 (D) -x 5 _ _ . A 
x 5 (D) ’°’ U2 W) 

uf (X) := (p , 1 > U 2 (X) 

where u 2 (x) is given by ( 21 ). 

• Compute uf D (x) according to: 


Ii^ CD (x) := < 


<P 


<P 


f v max v 
x 6 ~ x 6 

v max 


x 6~ x e 

I v min | 


, 0 ,u 5 (x)j if x 6 > x\ 
,l,^(x)j ifx 6 <x( 


( 22 ) 

(23) 


where u\ (x) is given by (17). 

• Compute the abscissa x 5 of the horizontal projection of p on BCD 

• Compute the vertical interpolation u£(x) according to 


« 2 (X) := { 


<P 

(P 


r v max 

X R ~ x 6 


X 


max 


x 6 ~ x c 


l ,min | 


o) 


if x 5 < x 5 
if x 5 > x 5 


• The feedback u 2 = I< 2 (x | P L ) is obtained by horizontal interpola¬ 
tion between u£(x) and uf CD (x) according to: 

^,uf D (x),u|(x)) (24) 


K 2 (X I P L ) := (p 


x 5 


5. The state feedback law 


Note that by definition of the interpolated feedback laws 
J<i(x,P L ) and J< 2 (x,Pl), the resulting closed-loop trajectories passes 
through neighborhoods V 34 (X 3 ,X 4 ) and V 5 6 (Xg,Xg) of the targeted 
pairs (X 3 ,x^) and (x^,x^) of size 0(1/2) after a finite time. In other 
words by choosing 2 > 0 (used in the interpolation map (8)) suffi¬ 
ciently high, the trajectories cross neighborhoods that are as small 
as required. In particular, 2 > 0 can be taken sufficiently high so as 
to make these neighborhoods entirely contained in the basins of 
attraction of the locally asymptotically stable dynamics under 
the local control laws uf(x^) and uj (x). If then, the feedbacks 
uf (X 3 ) and 1*2 (x) are maintained, then asymptotic stability follows. 

The preceding discussion suggests a dual-mode controller. 
Namely, when the state (x 3 ,x 4 ) [resp. (x 5 ,x 6 )] is far from the de¬ 
sired pair, the interpolated law J <3 [resp. I< 2 ] is used to steer the 
state to small neighborhood of the desired pair. When the state 
is inside this neighborhood, a switch to the locally asymptotically 
stabilizing law uf (x£) [resp. uj(x)] is operated. 

It is worth underlying also that if the neighborhoods are taken 
sufficiently small (that is 2 is taken sufficiently high), then the tran¬ 
sient excursions still meet all the saturation constraints on the 
states. The overall control law is therefore admissible. 

Note however that such a dual-mode control design has to be 
done with care in order to avoid undesirable chattering-like behav¬ 
ior (the state enters and leaves the small neighborhood with a high 
frequency). This is generally avoided by defining the size of the fi¬ 
nal neighborhoods using the Lyapunov functions level of the local 
feedback. This is more rigorously defined as follows: 


Definition 4. For each admissible X 3 , the Lyapunov matrix reflect¬ 
ing the stability of A(uf (x 3 )) is denoted by S 34 (x 5 ). Moreover, for 
any given pair (x 3 ,x 4 ), the following function is defined: 


d 34 (x,P L ) 


0 \ 
0 

x 3 -X3 

\x 4 -x s J{xl)J 


Fig. 5. Typical disposition of the vector fields of the dynamic system (le) and (If) 
for a given x 4 and for x 7 > x r 5 . 


S 34 (x*) 


(25) 
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where x* 3 is given by ( 20 ). 


Definition 5. Define S 56 e U 2 to be the Lyapunov matrix corre¬ 
sponding to the following Hurwitz matrix: 


(S? aw ) 

Vo -pj 



Moreover, for any pair (x 5 ,x 6 ), the following function can be 
defined: 


d56 (X,Pl) 


f X 5-^5 
\x 6 -x r 6 (x,P L ) 


$56 



where x r 6 is defined by (13). 

By now, we have everything we need to define the complete 
state feedback: 


Definition 6. Let 2 > 0 be the positive scalar invoked in the 
definition of the interpolation map (8). Let y 34 > 0 and y 56 > 0 be 
some positive reals. Let P L be a load power demand. The state 
feedback law is defined as follows: 


Ui := 


I<i(x,P L ) if d 34 (x,P L ) > y 34 

uf(Xj(x,P L )) otherwise 

(28) 

K 2 (x,P l ) if d 5S (x,P L ) > y 56 
u* 2 (x,P l ) otherwise 

(29) 


u 2 


Based on the discussions above, the following result can be 
proved: 


Proposition 1. If Assumptions 1-6 are satisfied, then there exist 
sufficiently small y 34 > 0 , y 56 > 0 and sufficiently high 2 x max 
( 734 , 756 ) > 0 such that the closed-system obtained by applying the 
state feedback laws (28), (28) to the dynamic system ( 2 a)-( 2 f) and 


( 2 g) asymptotically stabilizes the system’s state at a steady value that 
is compatible with the desired values x r 5 and x r 7 of x 5 and x 7 
respectively while meeting the problem’s constraints (3) and using 
admissible control u e [ 0 , l] 2 . 

Proof. See Appendix D. □ 

6. Simulation-based validation 

In this section some simulations are proposed to illustrate 
several facts regarding the proposed state feedback design. More 
precisely, the validation is split into two stages: 

1. In the first stage, referred hereafter as the decoupled simulations, 
each feedback law is validated separately by decoupling the two 
subsystems. This is done by simulating: 

• the system with state (xi,x 2 ,x 3 ,x 4 ) t in which x 5 is 
used, 

• the system with state (x 5 ,x 6 ,x 7 ) with a given and constant 
value of x 4 . 


These decoupled simulations enable us to show the impact of 
the saturation on the different state variables on the settling time 
of the corresponding controlled system to be shown. This stage 
obviously shows that there are two separate time scales which 
explains the success of the whole scheme. 

2. In the second stage, the whole closed-loop system is precisely 
simulated under varying power demand in order to assess the 
ability of the control law to tightly regulate the bus voltage x 5 
under varying power demand P L while regulating in a slower 
time scale the super-capacitor voltage x 7 . 
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Fig. 6. Control of the (x 5 ,x 6 ) system under constant x 4 . Influence of the allowable values of x 6 on the settling time. Comparison between the cases x™ ax = 50 and Xg ax = 60. 
Note that the time axis is in milliseconds. The closed-loop system shows a minimum-time like behavior since it saturate the constraints on the current x 6 = I L bb in order to 
accelerate the convergence. 
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6.1. Decoupled simulations 6.2. Global simulation 


Fig. 6 shows how the voltage x 5 is regulated around the de¬ 
sired value x d 5 = 50 V while the desired value x d A = 14.7 is consid¬ 
ered. Two scenarios are plotted in this Figure for two different 
values of the maximum current x™ ax e {50,60}. These scenarios 
clearly shows that closed-loop behavior shows a sort of mini- 
mum-time behavior since it saturates the current in order to 
accelerate the convergence towards the desired value. Fig. 7 
shows how the state x 4 (and therefore x 3 ) are steered to the de¬ 
sired steady pair which corresponds to + 20 % increase in 

the value of x d A while meeting the constraints on all the state var¬ 
iable. In particular, a minimum value of x% in = 4 A. This figures 
clearly shows how the bounding values defined by (A.2), (A.3) 
which guarantee the fulfillment of the constraints on x 2 slow 
down the possible evolution of the pair (x 3 ,x 4 ). Fig. 8 shows the 
same kind of results when a decrease of the desired value x d A of 
-20%. Note that in this case, the upper bound x% ax = 4.5 A is sat¬ 
urated in order to accelerate the convergence. It is worth noting 
that during the scenarios of Figs. 7 and 8 , the state variable x 4 it¬ 
self does hit the constraints during the first milliseconds. This is 
difficult to see because of the time scale. This is the reason why 
the same evolutions are plotted in Figs. 11 and 12 using a loga¬ 
rithmic time scale. 

Figs. 9 and 10, shows the impact of the saturation level x™ n and 
x™ x on the state variable x 2 on the response time of the closed-loop 
system when ±20% change in the desired value x d A is applied. Flere 
again, the closed-loop shows a sort of minimum time behavior 
since it systematically saturates the constraints in order to acceler¬ 
ate the convergence. 


x\ = Q (rad/sec) 



The decoupled set of simulations proposed in Section 6.1 clearly 
shows that the time scale of the two subsystems are quite sepa¬ 
rated. Indeed, the response time of the bus voltage related subsys¬ 
tem is in the range of milliseconds while the response time of the 
current x 4 = I Lfb that is used to charge the super-capacitor voltage 
is in the range of seconds. This justifies the decoupled design. Nev¬ 
ertheless, the global control design is validated in the present sec¬ 
tion through closed-loop simulation of the whole system under 
varying power demand scenario. 

Figs. 13-15 show three different scenarios that investigate the 
overall closed-loop performance under time varying power load 
demand profiles. In particular, Figs. 13 and 14 show two scenarios 
that differ only by the maximal allowed value x^ ax on the current 
x 2 = I red which is set to x™ ax = 5 A and x% ax = 4.8 A respectively. 
Note how this slight difference has a great impact on the amplitude 
of the excursion of the supercapacitor voltage x 7 = V sc . This sug¬ 
gests that the maximum allowable current x 2 which heavily 
depends on the choice of the component has to be jointly designed 
with the choice of the supercapacitor. 

All the simulations clearly show that the control objectives are 
met in the sense that the bus voltage is tightly regulated while the 
supercapacitor voltage is restored as soon as the constraints on the 
different control and state variables enable it. 

7. Conclusion 

In this paper, a global constrained control framework is pro¬ 
posed for the stand-alone hybrid (Stirling engine)/supercapacitor 
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Fig. 7. Control of the (x 3 ,x 4 ) system under constant x 5 = X 5 after a step change of +20% on the desired value x A . Note that the bounds x 3 and x 3 slow downs the response in 
order to meet the constraints on x 2 . The same scenario is depicted in Fig. 10 using logarithmic scale in order to show the evolution of the variables during first milliseconds of 
the scenario. 
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Fig. 8. Control of the (x 3 ,x 4 ) system under constant x 5 = Xg after a step change of -20% on the desired value x 4 . Note that the bounds x 3 and x 3 slow downs the response in 
order to meet the constraints on x 2 . The same scenario is depicted in Fig. 11 using logarithmic scale in order to show the evolution of the variables during first milliseconds of 
the scenario. 


power generation system. The control architecture explicitly han¬ 
dles saturations on the control variables as well as on the compo¬ 
nents of the state vector. The resulting closed-loop behavior can be 
used in the design step where the choice of the components and 
the class of power demand profiles have to be investigated. This 
is because such choices highly impact the size of the supercapaci¬ 
tor that can afford the different power demand profiles. 

Beyond the specific system that has been considered in the 
present paper, it is our claim that the control design methodology 
can be extended to a wide variety of power conditioning stages 
that can be encountered even when other thermodynamic cycles 
and/or different storing devices are used. 

Appendix A. Proof of Lemma 1 


modes and this regardless of the value of the constant steering sig¬ 
nal x 3 . This is because after this value t\ the derivative of y becomes 
monotonic and cannot change its sign anymore. 

Now for this once-for-all fixed value t*, the expression of the 
output y(t*) = x 2 (t*) can be given in terms of the current value of 
(xi,x 2 ) and the input x 3 using constant and off-line computable 
matrices M*,N* and IC according to: 

x 2 (t*) :=M*. f * 1 ) +N , a 2 +K’’x 3 (A.l) 

This obviously leads to the following expressions for x 3 and x 3 : 

j - N"a 2 


x 3 (xi,x 2 ) := max 


yinin 

^3 5 


X 


max 



Let us denote by z f and z s the fast and the slow eigen-modes of 
the linear system (2a), (2b) and denote by y = x 2 the output we are 
interested in. By definition of the eigen-modes, the following 
expression is derived: 

y(t) = oc f Zf(t) + oc s z s (t ) 

Suppose that we have some current initial state that satisfies the 
saturation constraints (3) on x 2 . Let us characterize all constant pro¬ 
files on x 3 leading to trajectories that still meet these constraints. 
According to Assumption 1 on 2/ and 2 S , the temporal generic struc¬ 
ture of the trajectory y( ) is approximately the one depicted in 
Fig. 16 (The first order behaviors of the eigen-modes are replaced 
by straight lines). Namely, the position in time of the singular point 
(zero time-derivative) is close to the settling time (t*) of the fastest 


x 3 (xi,x 2 ) := min 



v mm 

a 2 




(A.3) 


Now if x 3 is not constant but satisfies at each instant the constraint 
(4), the result can be proved by obvious receding-horizon way by 
looking at the first order two eigen-systems. This obviously ends 
the proof of Lemma 1. □ 


Appendix B. Mathematical conditions to check items(5), (6) and 

(7) of Assumption 4 

The Item (5) of Assumption 4 can be checked through the fol¬ 
lowing conditions that can be checked by scalar optimization in 
the decision variable x 3 : 
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Ired (A) 



^4 = ILfb (A) 



Fig. 9. Control of the (x 3 ,x 4 ) system under constant x 5 = Xg after a step change of 
system. 
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20% on the desired value x 4 . Influence of the constraints on x 2 on the response time of the 
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Fig. 10. Control of the (x 3 ,x 4 ) system under constant x 5 = Xg after a step change of +20% on the desired value x 4 . Influence of the constraints on x 2 on the response time of the 
system. 
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x\ = Q (rad/sec) 
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Fig. 11. Same scenario as the one depicted on Fig. 6 using logarithmic scale on the time axis in order to show the beginning of the scenario. This scale enables the saturation 
on x 4 to be clearly shown. 


min 

(X 3l X™)€AB 

max 

(x 3 ,x™ ax )eDE 


4 ( 4 ) 


4 ( 4 ) 


x,m\\\x,r 
A 4 a 5 

*3 


v max v r 
a 4 a 5 


*3 


>0 

(B.l) 

<0 

(B.2) 


since the analyzed expression is nothing but the resulting derivative 
of x 3 when x 2 , x 4 and x 5 take the appropriate values. 

As for the items (6) and (7) regarding the atractivity of the set 
BCD, it can be checked through the following conditions (in which 
the notation p = (x 3 ,x 4 ) r is used for compactness): 


min d(p(t,p°,Ui = 1), ABCDE) = 0 

(tjP)e[0,t f ]xn M v 7 

(B.3) 

min d(p(t,p°, iq = 0), ABCDE) = 0 

(t,p 0 )6[0,t f ]x£ 34 

(B.4) 


Step 1: We use the variable x 6 as a virtual control to regulate 
x 5 =V bU s around x r 5 . 

The regulation error being e 5 =x 5 -x r 5 , we use the following 
Lyapunov function: V 5 = \ ■ e\. Then: 

Vs = 65 • 65 = 65 • (aio • (X4 + *6) - — 1 — ■ -p -) 

'linv A 5 

By setting V 5 = -p 5 ■ ej < 0, with p 5 > 0 as a design parameter, we 
get: 



Pl 

Pond 


X 4 +^- (X 5 
a-io 



Step 2: To regulate x 6 around Xg e/ , the control variable u 2 is used. 
The regulation error being e 6 = x 6 - Xg e/ , we use the following 
Lyapunov function: V 6 = V 5 + \ ■ e\. Then: 


where the subsets 7 Z e R 2 and £ e M 2 are the sets that denote the 
values of (x 3 ,x 4 ) that lie to the right and to the left of BCD respec¬ 
tively. p(t,p°, Ui = 1) denotes the trajectories of (2c)-(2d) starting 
from p° under the control U\ = 1 and when x 5 = x^ and 
x 2 = (x^) are used. Similar definitions hold for p(t,p°,iii = 0) ex¬ 

cept that the control U\ = 0 is used. The notation d(p, ABCDE) de¬ 
notes the distance between p and the curves ABCDE. Finally, the 
final time t f is some sufficiently long computation time in order 
for the trajectory to cross the curve ABCDE. 


XT * • , , ypf , tl'JQ P/. 

Vq = 6s • Cs A- • 6 q = 6s • aio • X 4 + aio • (x 6 J + -) 

*lmv *5 

+ e 6 • (-On • x 5 + On • x-i ■ u 2 - x 6 ref ) 

= -p 5 • + +e 6 • (-a„ ■ x 5 + a„ • x 7 • u 2 - x 6 re/ + a 10 • e 5 ) 

By setting V 6 = -p 5 • 65 - p 6 • ej < 0, with p 6 as another design 
parameter, we get: 

«2(*) : = 7, “TT" ' ( Q H ' *5 - 0,0 ' (X 5 - X^) + xf - p 6 • (X 6 - Xf)) 

Qn • X 7 


Appendix C. Proof of Lemma 3 

To locally regulate x 5 =V bus around x r 5 , A backstepping design 
method is applied to Eqs. 2e,2 . 


Appendix D. Proof of Proposition 1 

Note that by definition of the interpolated feedback laws 
/Ci (x,P L ) and I< 2 (x,P L ), the resulting closed-loop trajectories passes 
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xi=Q (rad/sec) i 2 = Ired (A) 
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Fig. 12. Same scenario as the one depicted on Fig. 7 using logarithmic scale on the ti 
on x 4 to be clearly shown. 
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Fig. 13. Performance of the overall closed-loop system under varying power load demand. Case where 
in Fig. 13 where the upper bound x™™ = 4.8 is used. 


5. This scenario is worth comparing with the scenario depicted 
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X\ = fl (rad/sec) and Pi (W) X2 = Ired (A) 



%5 = Vbus and x 7 = V sc (V) 


u\ and U 2 



Fig. 14. Performance of the overall closed-loop system under varying power load demand. Case where x™ 1ax = 4.8. This scenario is worth comparing with the scenario depicted 
in Fig. 12 where the upper bound x™ ax = 5 is used. 
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X 5 = Vbus and x 7 = V sc (V) 
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Fig. 15. Performance of the overall closed-loop system under sinusoidal power load demand. 


through neighborhoods and V 56 (x 5 r,x^) of the targeted 

pairs (x^xj) and (x£,x£) of size 0(1/2) after a finite time. 

Let first y 34 > 0 and y 56 > 0 be taken sufficiently small to guar¬ 
antee that the corresponding level sets belong to the regions of 
attraction of the local controls uf(x5(x,P L )) and u^x.Pl). Once this 
is done, take 2 > 0 sufficiently high to guarantee that each of the 
conditions d 34 (x,P L ) < y 34 and d 56 (x,P L ) < y 56 will be reached in fi¬ 
nite time so that the local controllers can be fired. 

Now since the switching conditions in (28), (29) are defined 
using the Lyapunov level sets associated to the local stabilizing 


control laws, no chattering behavior occurs and the state is asymp¬ 
totically steered to the stationary state x st defined by: 



(D.l) 

4' = 4 

(D.2) 

4 = o 

(D.3) 

v st _ v r 

Ay - Ay 

(D.4) 

where x\ is given by: 
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Fig. 16. Illustration of the proof of Lemma 1: Under the conditions of the Lemma, 
the position of the pic of y = x 2 is almost independent of x 3 . 



Pl 

Pinv X 5 


(D.5) 


Indeed, this last equation comes from (14) because of the asymp¬ 
totic convergence of x 7 towards x r 7 that results from (15). This with 
(13) shows that x 6 converges to 0. This obviously ends the proof. 

□ 
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